Magnetic-field switchable metal-insulator transitions in a quasi-helical conductor 
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We study Anderson localization in disordered helical conductors that are obtained from one- 
dimensional conductors with spin-orbit interaction and a magnetic field, or from equivalent systems. 
We call such conductors "quasi-helical" because the spins of the counterpropagating modes are not 
perfectly antiparallel and have a small spin-wavefunction overlap that is tunable by the magnetic 
field. Due to the overlap, disorder backscattering is possible and allows a localization transition. 
A conductor can pass through two localization transitions with increasing field, one from the con- 
ventionally localized system to the quasi-helical conductor (with localization length exceeding the 
system length), and one at a higher field again to a localized state, due now, however, to backscat- 
tering below the magnetic field induced pseudo-gap. We investigate these transitions using a unified 
two-step renormalization group approach. 



I. INTRODUCTION 

Over 50 years ago, Anderson showed that a metal- 
insulator transition can arise due to localization of parti- 
cles by scattering on a disorder potential.— Since then, 
Anderson localization has evolved into an important 
topic of condensed matter physics, photonics, and ultra- 
cold atom gases.— While the basic localization mechanism 
can be understood on a single-particle basis, a complete 
investigation needs to include further interactions, espe- 
cially if they compete against localization. To study such 
a competition it would be advantageous to control the 
interactions externally, and to pass through the metal- 
insulator transition on demand. Here we show that this 
can indeed be achieved in the presently much investi- 
gated one-dimensional (ID) helical conductors through 
an external magnetic field. 

Helical conductors are characterized by spin-filtered 
transport, in which opposite spins (or Kramers part- 
ners) are bound to the right (R) and left (L) moving 
conduction modes. In conventional ID conductors, even 
weak disorder is already sufficient to turn the conductor 
into an insulator—^ by disorder-induced backscattering 
between the Fermi points ±kp. In a helical conductor, 
however, backscattering is only possible together with 
a spin-flip, and the conductor is insensitive to normal, 
spin-preserving disorder scattering.— 

Helical conductors appear at the edges of topological 
insulators^ or in quantum wires^ or nanotubes^ with 
strong spin-orbit interaction (SOI). They have attracted 
much attention recently as they allow for spin-filtering, 7 
Cooper pair splitting^ and, if in contact with a supercon- 
ductor, the realization of Majorana end states^ While 
localization cannot occur in a perfect helical conductor, 
many of the investigated conductors are imperfect and we 
call them "quasi-helical" : the spins moving in opposite 
directions are not perfectly antiparallel. They provide 
the handle to tune localization externally. Examples are 



semiconductor nanowires with strong SOI in the presence 
of a uniform magnetic field 7 or, without SOI, of a spiral 
magnetic field,— which can also take the form of a spiral 
Overhauser field due to ordered nuclear spins.— 

The SOI shifts the spin f, I bands by the wave vectors 
±qo (Fig. [U dashed lines). A magnetic field B x perpen- 
dicular to the f) i axis has now two effects. It lifts the 
degeneracy at k — by opening a pseudo-gap (Fig. [T] 
solid lines), and it breaks time-reversal symmetry. Tun- 
ing the Fermi level to the center of the pseudo-gap by let- 
ting kp = go allows conduction only through the modes 
close to momenta ±2go with opposite spins. Through B x , 
however, these spins are no longer antiparallel, disorder 
backscattering becomes again possible, and localization 
can occur. 

In this paper we provide a unified approach to such 
localization, taking into account disorder, magnetic field, 
spin overlaps, and electron interactions. We formulate 
a two-step renormalization group (RG) approach within 
the Luttinger liquid (LL) framework that provides us 
with a transparent picture of the underlying physics. We 
consider the zero temperature case, valid if the tempera- 
ture is smaller than any of the gap sizes. 

As a result, we find that varying B x can cause local- 
ization transitions of two kinds. A localized system 3 ' 4 
at B x = can make a transition to a quasi-helical con- 
ductor at a critical field B* lying (for the example of 
InAs nanowires) in the range up to ~ 1 T. This transi- 
tion is quite abrupt and appears when the B x generated 
pseudo-gap overcomes the disorder gap, and the quasi- 
helical state is lower in energy. The transition is marked 
by a rapid increase of the localization length £i oc , and 
the system is conducting if £i oc is greater than the sys- 
tem length L. Through the equal-spin overlap, £i oc then 
becomes strongly B x dependent and decreases with in- 
creasing field. At a larger critical field B x * > 1 T, when 
£ioc < C, backscattering between ±2g causes again lo- 
calization, resulting in a 4go modulated density wave co- 
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FIG. 1. Electron dispersion for the noninteracting quan- 
tum wire with SOI for B x = (dashed lines) and B x 7^ 
(solid lines). The magnetic field opens a pseudo-gap Az ~ 
fiBg\B x \/2 at k = 0. We denote the resulting lower (upper) 
band by a (6). At any energy E the left (L) and right (i?) 
moving states are labeled by \a/b, L/R) as indicated in the 
figure. The color gradient from red to blue indicates the spin 
mixing by B x . 



existing with the transverse spin polarization generated 
by B x . This coexistence of density wave and uniform 
polarization distinguishes this phase from conventional 
localization, and we shall call it "sub-gap localization." 
For very strong disorder, the helical phase is suppressed 
(B* — B**), and a transition takes place directly be- 
tween the two localized phases. We show concrete re- 
sults for InAs nanowires in Sec. IVII below. The under- 
lying physics, however, is general and applies to many 
materials, except for the edge states of topological in- 
sulators due to their different band structure. For the 
latter, quasi-helicity can still be obtained, but has differ- 
ent interesting consequences^ 

The plan of the remainder of the paper is as follows. 
In Sec. |TT] we introduce the model of the quasi-helical 
conductor. The approach for its solution in the presence 
of disorder is discussed in Sec. IIIII In Sec. IIVI we present 
the necessary background of the bosonization framework, 
and in Sec. IVlwe discuss the technical details of the two- 
step RG approach. Section IVII contains the discussion 
of the results for the example of InAs nanowires and the 
conclusions. 



II. MODEL 

We consider a generic interacting ID quantum wire 
of length C with SOI and an external magnetic field, 
described by the Hamiltonian 
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FIG. 2. Overlap integrals | (a, R \a, L) | 2 and | (a, R \b, L) | 2 
as functions of energy E > Az = ^Bg\B x \/2 w 0.5 meV 
(for InAs with B x = 2 T). 11 The states \a/b,L/R) are as 
indicated in Fig. [TJ The inset shows the B x dependence of 
(a, R\a,L) | 2 at energy E — fi in the middle of the Az gap. 



Here ^> s (x) are the electron operators at position x for 
spin s = t,| = +, — , p = —ihd x is the momentum oper- 
ator, (i the chemical potential, m the band mass, cxr the 
Rashba SOI strength, and Az = HBgB x /2 the Zeeman 
interaction strength, with Bohr magneton fi b and Landc 
g-factor g. U e - e is a general electron-electron interaction 
and Udis the disorder potential. a x ' z are the spin Pauli 
matrices, with the spin axes chosen such that an couples 
to a z and B x to a x . Equation (fTJ) describes a quantum 
wire with a single transverse subband, and we exclude 
the influence of higher subbands. 14,15 Realizations of Eq. 

© are found in GaAs, InA Sf ±i>i&±£ Ge/Si^^ or InSb 

20 

and Udis, Eq. (JTJ leads to the bands 
[[J The SOI shifts the t4 disper- 
ztmafl/fi.. The B x field opens a 
pseudo-gap at momentum k — by spin-flip scatter- 
ing. We denote the lower (upper) resulting band by a 
(b). Diagonalizing H leads to the dispersions E^' b = 
H 2 k 2 /2m, =p yj A 2 Z + h 2 a 2 R k 2 , and the wave functions 

\a,k) = u_ k |t) - u k ||), \b,k) = u k |t) + U- k It), with 
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the form (a/b, k\a/b, k') weight the disorder backscatter- 
ing amplitudes at a given energy. Figure [2] shows the 
two overlaps which will be relevant for the localization 
transition. 
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III. APPROACH 

To study a disordered interacting conductor in ID, it is 
generally convenient to bosonize Eq. (JT]) and use an RG 
approach. For the present model, however, this approach 
exhibits some peculiarities. First, the overlap integrals 
depend on k and the problem is no longer scale invari- 
ant, a central assumption for RG. Second, the low-energy 
physics below the scale set by Az is described mostly by 
the band a. Hence, the RG procedure of reducing the en- 
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ergy scale can lead to a problematic crossing the bottom 
of the b band, which also invalidates the bosonization 
formulation that neglects the band curvature. While the 
exact treatment of every aspect would require a sepa- 
rate investigation, the following two observations allow 
us to still implement a global RG scheme capturing the 
relevant physics in a transparent way. 

First, Fig. [2] shows that the spin overlaps above Az 
approach their asymptotics only slowly with increas- 
ing energy E (as | (a,R\a,L) | 2 - A%/2ma%E and 
| {a,R\b,L) | 2 ~ 1 - A 2 Z /4E 2 ). Considering constant 
overlaps evaluated, e.g., at half of the high energy cutoff 
scale yields a valid, rather conservative estimate of the 
true influence of the overlaps. Second, the interactions 
also renormalize Az and it grows quickly under the RG 
flow . 11 ' 12 Hence the reduced RG bandwidth E meets the 
growing Az at a value E — A* z 3> Az well above the 
bottom of the b band, at which the linearity of the bare 
dispersion remains valid. 

Due to these two properties, we propose a unified two- 
step RG approach to the effect of disorder based on the 
LL theory. In the first step, we integrate over high ener- 
gies far above Az using a constant backscattering overlap 
and the standard inclusion of the disorder potentiali 3 -^ If 
the disorder is strong enough, localization already occurs 
in this regime by the conventional backscattering mech- 
anism. Otherwise, at E = A* z , we proceed to the second 
step. 

Crossing the gap A* z corresponds to freezing out the 
interaction-generated density fluctuations that renormal- 
ize Az, and hence is not a singular transition through the 
true band bottom. While the description of the proper 
transition may be quite challenging, the physics below 
A* z becomes again simple. It is described by a differ- 
ent LL theory for the modes originating from the a band 
only 11 ' 12 ' 21 In this regime it is legitimate to use a larger 
spin overlap value and, within the accuracy of LL and 
RG theories, we choose to take the overlaps at the chem- 
ical potential /i (see Fig. [5J inset). The evaluation of 
the disorder backscattering (the second RG step) follows 
then the standard lines of a spinless LL, and localization 
occurs if £i oc < C. 



IV. BOSONIZATION 

Above the pseudo-gap, all fluctuations of the a and 
b bands must be taken into account. As we account 
for the overlap integrals separately, the standard ,s =f 
, i basis rather than the a, b basis 22 is more conve- 
nient for bosonizing. The electron operators are de- 
composed into R and L moving components as ip s (x) — 
e^'^fil+e-""' 1 ^,.^), for k FRs ,k FLs the two 
Fermi points of the (nominal) spin s bands. We then 
write ip r}S (x) = rj r ^ s exp (riT/n[(p s (x) + ri9 s (x)}) /V2nfi 
for r = R,L = +,— , the boson fields ip s and # s sat- 
isfying [(p s /(x'),d x '& s (x)] = i5 ss >S(x — x'), the Klein 
factors ?7 r . s , and the short distance cutoff k. Defin- 



ing <^ p ,i? CT = (<£t ± (fli)/V2 and <& p ,<p„ = (#t ± 
t9^)/V2, we obtain the Hamiltonian 22 (for Udis = 0) 

where K pa measure the interaction strengths (0 < K p < 
1 for repulsive C/ e -e and 22 K a « I) and u p i<7 are renor- 
malized velocities. Because /i lies in the pseudo-gap, Az 
renormalizes as well, expressed by a relevant Hamilto- 
nian oc Az cos(v / 27r(0p + &a)), acting on the fields of the 
b band . 11 ' 12 We neglect further existing Cooper scatter- 
ing processes^ which are overruled by the renormaliza- 
tion of Az and the disorder here but otherwise would 
dominate the physics. 

The second RG step starts when the growing Az 
meets the reduced bandwidth E at E = A* z . The 
6-band fields are then fully gapped ) 11 ' 12 and the low- 
energy theory is described by the fields related to the a 
band alone. The corresponding Hamiltonian is obtained 
from H bos by suppressing all 6-related fields, which leads 
Hg os = {u a /2)$dx[K- 1 {d x v a ) 2 + K a {d x V a ) 2 }, 
with u 2 a = ^/^[ul + ul+UpUeiKpKe + K^K- 1 )} 
and K 2 a = KpK-^UpKp + U<7 K- l )/{u p K- 1 + u a K p ). 
We have neglected here also a marginal coupling between 
a and b fields^ 2 - because it affects mostly the fermionic 
respons e 21 ' 23 but not the bosonic theory. 



V. DISORDER AND RENORMALIZATION 
GROUP APPROACH 

Disorder can be expressed by a random potential Udis 
with Gaussian distribution 3 -^ that scatters between the 
bands i,j — (a/b, L/R). The scattering amplitude is de- 
scribed by the dimensionless disorder strength Dy , pro- 
portional to the square of the strength of each individ- 
ual scattering potential, which in turn is proportional 
to |(«|j)| 2 , hence Dij — \(i\j)\ 4 D. Following the stan- 
dard replica disorder averaging approach^ we obtain 
the scaling equations (including also the amplitude y ex- 
pressing bulk backscattering, see Ref. y) 

d l Kp = -u p K 2 (2D ab + D aa )/Au ai (2) 

d t K a = -K 2 a (D ab + y 2 )/2 + D aa /A, (3) 

diy = (2 - 2K a )y - D ab , (4) 

diD aa = (3-K p - K- 1 )^, (5) 

dif) ab = (3-K p - K„ - y)D abl (6) 

dmp = -K p u 2 p (2D ab + b aa )/A U(7 , (7) 

diu a = ~u a K a D ab /2 - u a K- z D aa /A, (8) 

d l 5(l)=[2-(K p + K- 1 )/2]S, (9) 

with I the running RG scale, D aa — -D( a ,i?,),(a,L)> and 

b a b = -D(a,_R),(h,L)- We have neglected scattering be- 
tween (b, R) O (b, L), (a,R) o (b, R), (a,L) o (6, L) 
after verifying that it has no effect. We have also defined 
8(1) = A z (l)/E(l) with E(l) = hv F /n(l) the running 
effective bandwidth, for k(1) = ne l and Fermi velocity 
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FIG. 3. Localization length £i oc of an InAs nanowire as a 
function of B x for D — 0.01 and £ = 5 fim. Close to B x — 
the system is localized due to conventional backscattering. At 
B x ~ B x « 0.9 T the system crosses over to a quasi-helical 
conductor and at B x ~ B" w 2 T to a localized phase due 
to sub-gap backscattering. The transition at B x results from 
the competition of the renormalization of Az and D a i, and is 
thus expected to be quite abrupt and independent of C. The 
transition at B x * is a strongly C dependent cross-over. 



The latter equations express the competition be- 
tween disorder backscattering and the delocalizing effect 
by repulsive interactions^ and by B x induced spin-flip 
scattering. Localization occurs when the latter is not 
strong enough such that D a b ~ 1 (D aa remains small 
above the gap) before we reach E = A* z (6(1) ~ 1) or 
k(1) > C. Otherwise we switch to the second step, de- 
scribed by H^ os , with the parameters K a and u a obtained 
from the resulting K p>ai u P u of the first step. As argued 
above, the spin overlap weighting D aa can now (discon- 
tinuously) be replaced by (a, L |a, R) evaluated at E = /i. 
The RG equations are 



diK a = ~K 2 a D aa /2, 
dib aa = (3-2K a )D aa , 
diu a = -u a K a D aa /2, 



(10) 

(11) 
(12) 



describing, for effectively spinless fermions, the competi- 
tion between localization and delocalizing superconduct- 
ing fluctuations^ The latter overrule disorder localiza- 
tion above the critical attractive interaction strength 3 
K a > 3/2. However, K a < 1 for repulsive interactions, 
disorder scattering is relevant, and localization occurs if 
D aa ~ 1 is reached while k(1) < C; otherwise the (finite) 
system is a helical conductor. 



VI. RESULTS AND DISCUSSION 
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FIG. 4. Phase diagram for InAs nanowires obtained from 
curves as in Fig. for various D. At D > 0.023 the quasi- 
helical phase is suppressed and there is a direct crossover be- 
tween the localized phases. Not visible on the shown scale is 
the localization threshold for B x — 0, occurring at D w 0.003. 



and expresses a smaller effective bandwidth). The local- 
ization length is £i oc = min^- n(l*) / Dij(l*) with I = I* 
the scale at which the RG flow stops (in the first or the 
second step). If a Dij(l*) = 1 is reached before n(l) > C, 
the system is localized, £i oc < C. We consider disorder 
strengths about D ~ 0.01, leading at B x = to £i oc ~ 0.3 
fim. For a sample length of, e.g., C = 5 /xm, the sys- 
tem is well localized. At small fields, the D ao and B x 
scattering processes compete, and if B x passes a critical 
value B* < 1 T, Az(l) overrules the disorder backscat- 
tering and the system becomes a quasi-helical conductor, 
where £i oc > C is now determined by the sub-gap disor- 
der strength D aa . At B x = B** ~ 2 T, \{a,R\a,L)\ 
becomes large enough such that & oc < C and the system 
crosses over into the sub-gap localized phase. In Fig. 
[3] we plot this crossover behavior for D = 0.01. Trac- 
ing similar curves for various disorder strengths D leads 
to the phase diagram shown in Fig. [4] For strong disor- 
der, the quasi-helical conduction phase is absent, and the 
crossover takes place directly between the two localized 
phases. 

To conclude, we note that the disordered quasi-helical 
system shows at B* a phase transition from the conven- 
tional Anderson localized phase in which £i oc only weakly 
increases with B x , to a phase in which £i oc jumps to a 
very large value and then decreases with increasing B x . 
For fields B* < B x < B**, £i oc exceeds C and forms 
a quasi-helical conductor before crossing over into an- 
other localized phase. While the former localized phase 
is characterized by the conventional 2kp = 2qo modu- 
lated charge density wave, the latter sub-gap localized 
phase combines a 4qo density wave with a uniform elec- 
tron polarization as the order parameter. The existence 
of a quasi-helical conduction state between the two insu- 
lating phases as a function of B x is a peculiar behavior 
that could be used to test if a conductor is helical. Re- 
markably the conducting phase extends into the regime 
of quite strong disorder, which indicates that a quasi- 
helical conductor does not necessarily require ultraclean 
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samples. The boundaries of this phase in the phase di- 
agram are controlled predominantly by A^, and so the 
g-factor of the material has the largest influence on the 
phase diagram. For instance, a similar phase diagram as 
Fig. g] for InSb with^S g m 50 spans over B, = 0-0.1T 
and D = - 0.002. 
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